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Abstract 



With (•) denoting an average with respect to the eigenvalue PDF for the Laguerre 
unitary ensemble, the object of our study is En{I', a, fJ-) (Yii'-iXfo oo)\ii^ ~ -^0^) 
for I = (0, s) and / = {s,oo), where X/' ~ 1 for Xi G I and Xj' ~ otherwise. 
Using Okamoto's development of the theory of the Painleve V equation, it is shown 
that _Ejv(/;a,/x) is a r-function associated with the Hamiltonian therein, and so can 
be characterised as the solution of a certain second order second degree differential 
equation, or in terms of the solution of certain difference equations. The cases ^ = 
and jj, = 2 are of particular interest as they correspond to the cumulative distribution 
and density function respectively for the smallest and largest eigenvalue. In the case 
/ = (s, oo), En{I', a, tj) is simply related to an average in the Jacobi unitary ensemble, 
and this in turn is simply related to certain averages over the orthogonal group, the 
unitary symplectic group and the circular unitary ensemble. The latter integrals are 
of interest for their combinatorial content. Also considered are the hard edge and soft 
edge scaled limits of En{I\ a, /i). In particular, in the hard edge scaled limit it is shown 
that the limiting quantity E^^^'^{{0, s); a, /i) can be evaluated as a r-function associated 
with the Hamiltonian in Okamoto's theory of the Painleve III equation. 
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1 Introduction and summary 

In a previous paper |17] the quantities 
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N 



1=1 



E^{X;a) :=(nX(-oo,A](A-AO'^ 
and 



GUE 



A(-oo,A] 



1 Xi e (-00, A], 
otherwise 



N 



:i.2) 



1=1 



GUE 
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where the averages are with respect to the joint eigenvalue distribution 
of the Gaussian unitary ensemble (GUE), were shown to be equal to 



the T-functions occurring in Okamoto's theory |29| of the Painleve IV 
equation. It was noted in |jl^ that we expect the analogous quantities for 
the Laguerre unitary ensemble (LUE) to be expressible in terms of the r- 
functions occurring in Okamoto's theory |3^ of the Painleve V equation. 
It is the purpose of this article to verify this statement by giving the details 
of the correspondences between the multi-dimensional integrals defining 
the analogues of E]\f{X; a) and Fn{X; a) for the LUE, and the r-functions 
from m. 

Let us first recall the definition of the LUE. Let X he a n x N {n > 
N) Gaussian random matrix of complex elements Zjk, with each element 
independent and distributed according to the Gaussian density ie"!^^*!^ 
so that the joint density of X is proportional to 

(1.3) e-T=^^'^. 

Denote by A the non-negative matrix X'^X. Because (|1.3| ) is unchanged 
by the replacement X i-^ UXV for U a. n x n unitary matrix and V a 
NxN unitary matrix, the ensemble of matrices A is said to have a unitary 
symmetry. The probability density function (PDF) for the eigenvalues of 
A is given by 

1 ^ 

(1.4) T^n^^""' n (Afc-A,f, A,>o, 

i=l l<j<k<N 

where C denotes the normalization and a = n — N,n>N. (Through- 
out, unless otherwise stated, the symbol C will be used to denote some 
constant i.e. a quantity independent of the primary variables of the equa- 
tion.) Because A'^e"'*' is the weight function occurring in the theory of the 
Laguerre orthogonal polynomials, and the ensemble of matrices A has 
the aforementioned unitary symmetry, the eigenvalue PDF (|1.4| ) is said 
to define the LUE. 

Analogous to ^Af(A;a) specified by 1^) for the GUE, we introduce 
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the quantities 

N 

(1.5) ^^((0, .); a, /.) := ( n x£^) (A. - ^T] 

N 

(1.6) En{{s, oo); a, /i) := ( I] xfl^{s - XiT] 

where the averages are with respect to ( |1.4[) (the parameter a in ( |l.lD has 
been denoted ^ in (^), (0) to avoid confusion with the parameter a in 
(|L4|)). Exphcitly 



(1.7) ^^((0,s);a,//) 



iV 

l<j<k<N 
Ns N 

n / dXi Af (A, + s)''e-^' n i^k - Xj] 



1=1-'^ l<j<k<N 



p -Ns iV oo 

>+M)A^+iv^ Jl / dA, Af (A, + l)"e-^^' n (A.-A,)' 

1=1-^^ l<j<k<N 



^ i=l-- l<j< 



(1.8) £'Ar((s,oo);o,/i) 

AT 



i n r ^^"'n^ - AO^ n (Afc-A,)2 

Z=l'^° l<j<k<N 
N 1 

n / dA,A?(l-AO^e-^^' n (Afc-Ai)^ 

/=l-^'^ l<j<k<N 



- ]_Ja+^J.)N+N2 

In the case ^ = the first integrals in ( [1.7] ) and ( |1.8D are the definitions 
of the probabihty that there are no eigenvalues in the intervals (0, s) and 
(s, oo) respectively of the Laguerre unitary ensemble. The case ^ = 2 also 
has significance in this context. To see this, first note from the definitions 
(using the first integral representation in each case) that 

(1.9) ^En+iHO, s); a, 0) a s"e-^^^((0, s); a, 2) 
as 

(1.10) ^^iv+i((s,oo);a,0) oc s''e~'EN{{s,oo);a,2). 
as 
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On the other hand, with Pmin('S;a) and Pmax(s;o) denoting the distribu- 
tion of the smallest and largest eigenvalue respectively in the N x N LUE, 
we have 

(1.11) Pmin(s; a) = --^^^((0, s); a, 0) 

as 

(1.12) pmax(s;a) = 3--^Af((s,oo);a,0). 

as 

Under the replacement N ^ N + 1, Pramis;a) and Praa.xis;a) are deter- 
mined by Eiy{{0, s); a, 2) and EnHs, oo); a, 2) respectively. 



From the second formula in (1.7), we see that with 

N 

/LUE 



(1.13) F^(s;a,/i):=(n(Ai-s^^\ 



1=1 



we have 

(1.14) Fn{s; a, //) = (e^^^7v((0, s); //, a)) 

(notice the dual role played by jj and a on the different sides of ( 1.14| )) 



so there is no need to consider F]sr separately. Note that with /x = 2, 



(1.13) multiplied by s^e"* is proportional to the definition of the density 
in the LUE with N -|- 1. We also remark that for /i non- integer 

the definition ( |1.13| ) must be complemented by a definite choice of branch 
of the function (A/ — s)^. But this case does not appear in our random 
matrix applications so we will not address the issue further. 



Some quantities generalizing (II^) and (|1.§| ), and which include ( |1.13 ) 



are 



(1.15) EN{{0,s);a,i2;0 

1 ^ P OO P s 

-■=rll{ -U )dXiXte~>^^{Xi-sr n i^k-Xjf 

l^-^ Jo Jo l<j<k<N 

(1.16) En{{s, oo); a, 

■■=rlli -d )dXiXte-''{s-Xir n i^k-Xjf. 

'^1=1^-^0 Js l<j<k<N 
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Thus 

£^Ar((0, s); a, fi; 0) = Fn{s; a, /i) 

En{{0, s); a, fi; 1) = £^{{0, s); a, /i) 

^Af((s,oo);a,/x;l) = EnHs, oo); a, fi). 

Note that only one of the quantities (|1.15| ), (|1.16|) is independent since 



EnHs, oo); a, fj.;C) 



V(l-A)^ 

Their interest stems from the facts that 



A=?/(«-l) 



71 ! 5^ 

(-1)" 9" 



-^7v((0,s);a,0,0 
^Ar((s,oo);a,0,0 



^7v(n; (0,s)) 
£'Ar(n; (s,cx))) 



n! 9^" 

where Ej\f{n;I) denotes the probabihty that in the LUE the interval / 
contains precisely n eigenvalues. 

The second integral in ( |1.8| ) is of interest for its relevance to the Jacobi 
unitary ensemble, which is specified by the eigenvalue PDF 

1 ^ 

(1-17) ^U^fi^-^if n (Afc-A,)^ 0<Az<l. 

1=1 l<j<k<N 

This ensemble is realised by matrices of the form A{A + B)~^, where 
A = X^X, B = Y^'Y, for X (Y) an m x N (n2 X N) complex Gaussian 
random matrix with joint density ( |1.3| ). The parameters a and b are then 
specified by a = rii — N , b = n2 — N (c.f. the value of a in (|1.4|) ). We see 
from ( |1.17D that 

(1.18) 



-n/ dxixfii-xi)'e^'' n 

«=1'^° l<j<k<N 



JUE 



so substituting in (|1.8D and compensating for the different normalizations 
in (jTI) and ( |lT7|) shows 

(1.19) ^^((.,oo);a,^) = M^.C^+Mm^Ve-Ef.! a, 

lN{a) \ 



JUE 
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where 



POO 

(1.20) /^(a):=n/ dXiXle-^' \{ {\k - \jf 

TV 1 

(1.21) jN{a, := n / - n - x,?. 

Z=l'^° l<j<k<N 

As is well known, the integrals /at (a) and Jn{cl, A*) can be evaluated in the 
form A^! Il^Jio^ '-i' where cj is the normalization of the monic orthogonal 
polynomial of degree j associated with the weight functions X"'e~^ and 
A'^(l — A)'' respectively. 

One important feature of the JUE is that with the change of variables 

(1.22) \j = ^{cos0j + 1), 

the PDF ( |1.17| ) assumes a trigonometric form, which for appropriate 
(A^, a, b) coincides with the PDF for the independent eigenvalues e*^^ of 
random orthogonal and random unitary symplectic matrices. In the case 
of orthogonal matrices, one must distinguish the two classes O'^ and 
according to the determinant equalling +1 or —1 respectively. The cases 
of N even and N odd must also be distinguished. For N odd all but 
one eigenvalue come in complex conjugate pairs e^*^^ , with the remain- 
ing eigenvalue equalling +1 for and —1 for Oj^. For N even, all 
eigenvalues of come in complex conjugate pairs, while for all but 
two eigenvalues come in complex conjugate pairs, with the remaining two 
equalling ±1. Let us replace A'^ in ( [1.17 ) by N* and make the change of 



variables ( |1.22 ). Then the PDF for the independent eigenvalues of an en- 



semble of random orthogonal matrices is (see e.g. |l^) of the form ( 1.17 ) 
with 



(1.23) 

(A^/2,-1/2,-1/2) for matrices in Oj^, A^ even 

((A^- l)/2, -1/2, 1/2) for matrices in O^, A^odd 

{{N - l)/2, 1/2,-1/2) for matrices in O;^, A^odd 

(A^/2 - 1,1/2,1/2) for matrices in O;^, iVeven 



{N*,a,b) 
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Matrices in the group USp{N) are equivalent to 2N x 2N unitary 
matrices in which each 2x2 block has a real quaternion structure. The 
eigenvalues come in complex conjugate pairs e^*^^ , and the PDF of the 
independent elements is of the form ( |1.17| ) with the change of variables 



L2^ ) and 

(1.24) iN*,a, b) = {N, 1/2, 1/2). 
It follows from the above revision that 

(1.25) (e^Tr{c/)\ 

for G = O^, or USp{N) is a special case of the more general average 
( p8|) . Exphcitly, 

,sTtiU) \ ^ -2sN* ^xn* Vp^'* S^f=i 



(1.26) e^^-'^^n =e-'''' e^^*'{e^ 

\ / U€G \ /JUE 

where on the RHS the dimension of the JUE is at first N* , then the 



parameters {N*,a,b) are specified as in (|1.23D or (|1.24| ). Also, xa^* = 
for G = USp{N) and O^, (A^ even), while xn* = ±1 for G = Ojj, 
(A^ odd) respectively. The averages ( |1.25| ) have independent importance 
due to their occurrence as generating functions for certain combinatorial 
problems. Let us quote one example. Set /^™^^ equal to the number of 
fixed point free involutions of {1, 2, ... , 2n} constrained so that the length 
of the maximum decreasing subsequence is less than or equal to 21, and 
introduce the generating function 

oo ,2n f(inv) 



[1.27) m)--=e-''''Y.l^ 



^„ 2"ra! (2n - 1)!! 

?i=0 ^ ' 

According to a result of Rains p3] (see also [S] ) we have 



t^/2/tTT{U)\ 



(1.28) m) = e ,,,,,.n. 

\ / U€USp{l) 

A list of similar results is summarised in ||2|. 

The second integral in ( |1.8D can be written in a trigonometric form 
distinct from that which results from the substitution ( |1.22| ). For this 
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purpose we recall the integral identity ^ 

N 



iV 1 

(1.29) n I dtitr' f{h,...,tN) 

. AT ^ /•1/2 
Vsinvre/ ,^J-i/2 



valid for / a Laurent polynomial. Applying ( |1.29| ) to ( 1.21 ) gives the 
trigonometric integral 

^ /•1/2 

(1.30) Mjv(a',6') := n / dx/ e^*^'(°'-''')|l + e2^^^'|'^'+^' 

1=1 "'-1/2 

X |g27rjxfc _ g27rjXj |2 

Moreover, applying ( |1.29| ) to the LHS of ( |1.18| ) and using Carlson's the- 
orem it follows that 



(1.31) (e 



JUE bi-^fj. 

N 

mxi{a'~b')\', I 2mxi\a'+b' se^'^^'^l 



Mn{0,0) 

MN{a',b')\f-J.'^ ' " ' /cue' 



J|g.*x,{a'-6')|i + g2 
1=1 



where a' = N + a -\- ^, b' = —{N + a) and the CUE (circular unitary 
ensemble) average is over the eigenvalue PDF 

(1.32) 1 Yl |e2™'= - e2™^ |2, -1/2 < x/ < 1/2. 

i<j<k<N 



Recalling ( 1.19| ) then gives 

N 

(1.33) / J|e™'('''-*')|1 + e2™'P'+*'e 



AT 

27ria:,|a'+fe'^-se2'^'^i 

1=1 'cue 



Mat (g , b )/jv(a) ^„ ^4-^^ jv_jv2 - 

MAr(0,0)JAr(o, /i) 

although the RHS is convergent for only a subset of the parameter values 
for which the LHS converges (explicitly the RHS requires Re(a' + 6') > — 1 
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and Re(— A^— 6') > —1 to be properly defined, while the LHS only requires 
Re(o' + b') > —1). Note however that by replacing the integrals over 



A; G [0, 1] in (|l.8|) by the Barnes double loop integral ||41[] the quantity 
E]y{{s,oo);a, fj,) has meaning for general values of a and fi. In the case 
a' = 0,b' = k£ Z>o 

N 

(1.34) (n(i+ 



the LHS of ( p:33D reads 

det(l + i7)^e" 



g— 27rjx; \fcg— se^'^" 



1=1 



CUE 



-sTi(U)\ 

I (7eCUE' 



Like ( 1.2^ ) this latter average with the dimension replaced by I is 
the generating function for certain combinatorial objects [^] (random 
words from an alphabet of k letters with maximum increasing subsequence 
length constrained to be less than or equal to /). 

As we have already stated, our interest is in the evaluation of (f]B) and 



1.7 ) as r- functions occurring in Okamoto's theory of Py • In some cases 



such evaluations, or closely related evaluations, are already available in 
the literature. One such case is (1.8) and (1.7) with fj. = 0, for which 
Tracy and Widom |Q (for subsequent derivations see [||, ^) have 
shown 



:i.35) 
:i.36) 



^7v((0,s);a,0) 
En{{s, oo); a, 0) 



exp 



exp 



a{t) 



dt 

a{t)_ 
t 



dt 



where a{t) satisfies the Jimbo-Miwa-Okamoto a- form of the Painleve V 
equation, 



;i.37) {ta 



with 



.//^2 



a-ta' + 2{a'f + (z^o + 1^1 + 1^2 + ^sW] 

+ 4(z^o + + cf'){v2 + fj'){v2. + = 0, 



;i.38) 



z^o = 0, 1^1 = 0, U2 = N + a, 



N 



(or any permutation thereof, since ( 1.37] ) is symmetrical in {fk}), and 
subject to appropriate boundary conditions in each case. For general ^ 
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the second integral in ( |1.8D has been written in an analogous form to 
the right hand side of ( |1.35D by Adler and van Moerbeke P], with the 
function corresponding to cr presented as the solution of a certain third 
order differential equation. It is remarked in Q that this equation can be 
reduced to a second order second degree equation, known to be equivalent 
to ( 1.37 ), using results due to Cosgrove |8| (see also the discussion in 
the final section of |42|). The average ( |l.34| ) can be computed in terms of 
Painleve transcendents by making use of an identity [^] implying 

(1.39) (det(l + U)''e-'^''-^A = e''' Ek((0, s); N,0), 

(note the dual role played by k and on the different sides of ( |1.3S| )) 
then appealing to ( 1.35| ) (direct evaluations have been given by Tracy and 



Widom 1 38 1 and Adler and van Moerbeke ||2[). 

The Okamoto r-function theory of Py provides a unification and ex- 
tension of these results. Consider first ( |1.7[ ). Proposition 3.2 below gives 



En[{^, s); a, /i) = exp / dt 

1n{cl) Jo t 



(1.40) 



In {a) 



IVa — N"? 

s e exp 



VN{t;a,^j.)+Nt-N{a-n) 



t 



dt 



where VN{t;a,^) satisfies the Jimbo-Miwa-Okamoto a form of Py ( 1.37 ) 
with 



(1.41) fo = 0, 1^1 = -fi, 1^2 = N + a, U3 = N, 
and subject to the boundary condition 

(1.42) y^(i;a,/x) ~ -Nt + N{a - ^l) - + Oih. 

t^oo t 

Regarding (1.8), in Proposition 3.1 below we show 



(1.43) 



In {a) 
s^^exp 



exp 



dt 



( 



UN(t;a,ii) 



dt 
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where like VN{t;a, UN{t;a,^) satisfies ( |1.37| ) with parameters ( [L.41 ), 
but with the boundary condition 



:i.44) 



t->0 a + IJ, + 21\ 



As well as VAr(t;a, /i) and Uj\[{t; a, fj,) satisfying the differential equation 
( |1.37| ), we show in Proposition |3.5| that they also satisfy third order dif- 
ference equations in both the a and fi variables. 

By substituting the first equation of ( 1.43|) in (|1.19|) and ( 1.33| ) we 
obtain the evaluations 



:i.45) 



JUE 



6i— >/x 



exp 



UN{t;a,ix) -aN - N"^ 



dt, 



N 



1=1 



CUE 



MN{a',b') [^UN{t]a,^l)-aN -N^ . 

exp / ^ dt 



M^(0,0) "Jo 

Comparison of ([L45|) with ([L2^ ) gives that for G = 0%, USp{N), 



:i.47) 



gXiV*s/e-sTr{(7) 



4s 



r*2 



= e--exp/ UM^-^-^b)-aN*-N* 
ugg Jo t 



with {N*,a,b) specified by ( 1.23D or (|l.24| ) as appropriate. In particular, 
with G = USp{N) and thus {N*,a,b) given by ( |1.24 ), substitution into 
(11381) shows 



^ .2/2 25/ f^'Ui{t;l/2,l/2)-l/2-f^ 

Pi(s) = e'" /^e^^' exp / ^ i ^ i > 1 

Jo t 



(1.48) 



After recalling the statement below (1.12) and inserting the propor- 
tionality constants in ( [L.9| ) and ( |1.10 ), the evaluations ( 1.40| ) and ( |1.43 ) 
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give 



(1.49) 

Pmin(^'5 

(1.50) 

Pmax(s;a) 



/Ar , .^ lN{a+2) , VN{t;a,2) + 2N 

= (A'+l)— -—s e exp / at 

N^N+i lN+i[a) Jo t 



= (iV+l)y^ 
N^N+l lN+i{a) 



(■ 



UN{t;a,2) 



dty 



Also, recalling ( 1.14 ) and the sentence below that equation, we see from 
( |OoD that 



(1.51) p{s) 

= iN+ 



^^7jv(a+2)^ag-(7V+i)5 
In+1 (a) 



exp 



V7v(t;2;a) + aN 



dt 



where p{s) denotes the eigenvalue density in the LUE. 



According to Proposition ^ 
(1.52) 



EmUO, s); a, 0; = exp / YMl^ dt 

t 



where W7v(t;a, 0) satisfies ( |1.37|) with parameters ( 1.41 ), fi set equal to 
zero, and boundary condition 



(1.53) 



WN{t; a, 0) ~ 



T{N + a + l) , 

^r(iv)r(a + i)r(a + 2) 



This result is known from |35]. 

A generalization of the integral identity ( 1.39| ) is derived in the context 
of the Painleve theory in Proposition below. We remark that this 
generalization is in fact a special case of a still more general identity, 
known from an earlier study |jll[|, involving an arbitrary parameter /3. 
For /? = 1 an average in the Laguerre orthogonal ensemble is related 
to an average in the circular symplectic ensemble, while for /3 = 4 an 
average in the Laguerre symplectic ensemble is related to an average in 
the circular orthogonal ensemble. 
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The Laguerre ensemble permits four scaled, large N limits. These are 

[0, i n, n 



(1.54) 
(1.55) 
(1.56) 
(1.57) 



a 
a 



s ^ s/4N, N ^ oo 

(7 - 1)A^, N{1- y/^f - u^{N)s, N ^00, 
(7 - 1)A^, s^N{l + + u+{N)s, N ^00, 



where in ( 1.56| ), ( |1.57D we require 7 > 1 and 
(1.58) 



z._(iV) 



iV(V7-l)(l-^) 



The first gives the limiting distributions at the hard edge, so called be- 
cause the eigenvalue density is strictly zero on one side. The second, third 
and fourth give the limiting distribution at the soft edge, so called because 
there is a non-zero density for all values of the new coordinates, but with 
a fast decrease on one side. 

The soft edge limit for the quantities ( |1.1D and ( |1.2| ) in the GUE was 
studied in our work [^], where evaluations in terms of solutions of the 
general Jimbo-Miwa-Okamoto cr-form of Pn 



0. 



(1.59) {u"f + 4u' ((n')^ - tu' + uj-a^ 

were given. From either of the limiting procedures ( 1.55D , ( |1.56 ) or ( |1.57 ) 
we reclaim the results of |l^ (because no new functional forms appear we 
will not discuss this case further in subsequent sections). In particular, it 
follows from the second equation in (1.43) that 



;i.6o) 



where 
(1.61) 

u{t,fi) 



lim Ce~'''/^EN{{s,oo);a,fi) 

sM4JV+2{2JV)l/3a 
JV-»oo 



^^°^*(so;/u)exp / u{t;fi)dt 



so 



lim 



AN [UN{t;a,^x) + Nii-^t 
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a ^ a — N ^ + —t, 



in ( p..37D with parameters ( 1.41[ ), and equating terms of order N"^ in the 
equation (which is the leading order) shows u{t; /i) satisfies ( |1.59 ) with 
a = /U. We know from [0] that ( 1.59 ) is to be solved subject to the 
boundary condition 



1/3, 



^2/3, 



(1.62) 



/ , 1 2 V _ 1 (4^2 - l)(4u2 - 9) 



As an application, ( 1.61 ) can be used to compute the scaled limit 
Pis) := limP,(/-i(20^/'s). 

(— >oo V Z ^ 

Imposing the condition that P{s) ^ 1 as s ^ cxd (which fixes the con- 
stants), a short calculation using ( |1.61 ) in ( 1.48| ) shows 

(1.63) P{s) = exp ( - J^{u{-t; 1/2) - t^/A) dt) . 

This is interesting because it has been proved by Baik and Rains Q that 



(1.64) 



P{2 



1/3. 



where Fi(s) is the cumulative distribution function for the largest eigen- 
value in the scaled, infinite Gaussian orthogonal ensemble. Tracy and 
Widom (see for a simplified derivation) have shown that 



(1.65) 



Fi{s) 



where with Ai(t) denoting the Airy function, q{t) is the solution of the 
non-linear equation 

(1.66) q" = 2q^ + tq, 

(Painleve II equation with a = 0) subject to the boundary condition 

(1.67) q{t) ~ -Ai(t) as t ^ 00. 
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By equating ( 1.63 ) and ( |1.64 ) we obtain as an alternative to ( 1.65| ) the 
formula 



:i.68) 



Fi (2-i/3s) = exp ( - y {u{-t; 1/2) - f/i) dt) . 



However the boundary condition ( |1.62| ) gives u (-t; 1/2) -tV4~0 which 
is vacuous, whereas the boundary condition according to the requirement 



Ai'(t) + (Ai(t))' 



The significance of ( 1.68 ) from the viewpoint of gap probabilities in ran- 
dom matrix ensembles as r-functions for Hamiltonians associated with 
Painleve systems will be discussed in a separate publication [^] . 

For the hard edge scaling ( [1.54 ), define the scaled version of (|1.7D by 



(1.69) 



y ( In (a) 
hm -— - 

s^s/4N \I]^[a + fl) 



In the case /i = 0, it is known from the work of Tracy and Widom |36] 
that 



(1.70) 



^'^^'^^(s;a,0) =exp 



t 



dt 



where (TB{t) satisfies the Jimbo-Miwa-Okamoto cr-form of the Painleve III 
equation 



1 



(1.71) {tcj"Y - viV2{a'Y + a'{Aa' - 1){(J - ta') - -^{vi - V2) 







with 



(1.72) 



vi = a, 



V2 = a. 



It follows from Proposition below that for general fi 

o-(t) + /i(/i + a)/2 



:i.73) ^^-'i(s;a,/x)=exp(-|J 



dt 
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where a satisfies ( |1.71 ) with 

(1.74) f 1 = a + ^, V2 = a — ^ 

and subject to the boundary condition 

t at'/^ " ' 



(1.75) 



+ 



V 4 2 



In Proposition 4^ we show that v{t;a,fi) = —cr{t) — fj,{iJ, + a)/2 satisfies 
third order difference equations in both a and fi. 
One application of ( |1.73 ) is to the evaluation of 



(1.76) 



hard 



;«)• 



First we note that the explicit evaluation of (|l.2[l|) is YljLi r(l+j)r(a+j) 
which together with the asymptotic formula r(x + a)/T{x) ~ x*^ for x — 
oo shows that 



hm ^ANW-f^s^e 



s^s/4N 22«+2r(a + l)r(a + 2) ' 
Using this formula and ( 1.73| ) we can take the large limit in ( |1.49| ) as 



required by (|1.76|) to conclude 



(1.77) <in(^;a) = ,..+.r(a + l)r(a + 2) 
where a{t) satisfies ( 1.71| ) with 

vi = a + 2, V2 = a — 2 



a(t) + a + 2 



dt 



and is subject to the boundary condition ( 1.75| ) with /x = 2. Similarly, it 
follows that the hard edge density is given by 



(1.78) p^^"^is 



exp 



a{t) +a{a + 2)/2 



22»+2r(a + l)r(a + 2) 
where a{t) satisfies ( |1.71| ) with 

vi = 2 + a, V2 = 2 — a 
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and is subject to the boundary condition ( 1.75| ) with (a,/x) = (2, a 



The result ( [L.73 ) also has relevance to the CUE. This occurs via the 



scaled limit of (1.39), which has been shown [^] to imply for a G Z>o, 



/i G Z, the identity 

(1.79) E^''''^{t;a,^i) oc e-^/H->"'l'^ie^^^'^'^^+^\detU)->') 



I c/ecuEa 



(eq. ( [4.32 ) below; the proportionality constant is specified in ( [4.33| )). It 



thus follows that this CUEa average can be evaluated in terms of the 
transcendent a{t) in ( 1.73]) (in the case /i = this has been shown directly 



in the previous works |3^, §). 

This concludes the summary of our results. In the next section (Sec- 
tion 2) an overview of the Okamoto r-function theory of the Painleve 
V equation is given. In particular a sequence of r- functions T[n], each 
of which can be characterised as the solution of a certain second order 
second degree equation, is shown to satisfy a Toda lattice equation. The 
fact that a special choice of initial parameters in the sequence permits 
the evaluations r[0] = 1 and r[l] an explicit function of t (a confluent 
hypergeometric function) then implies that the general member of the 
sequence r[n] is given explicitly as a Wronskian type determinant de- 
pending on t[1]. In Section 3 the Wronskian type determinants from 
Section 2 are evaluated in terms of the multiple integrals ( |1.7D and ( |1.8D , 
thus identifying them as r-functions and implying their characterization 
as the solution of a non-linear second order second degree equation. 

In Section 4 the Okamoto r-function theory of the Painleve III equa- 
tion is revised, and again a r-function sequence is identified which can 
be shown to coincide with a quantity in the LITE (explicitly E^'^'^'^{s\ a, fj.) 
for general fi and a E Z>o). By scaling results from Section 3 it is shown 
that for general /i and general a, E^^'^'^(s; a, /x) as specified in ( 1.69| ) is a r- 



function in the Pni theory. Concluding remarks relating to the boundary 
conditions ( 1.75| ) and (1.62) are given in Section 5. 



2 Overview of the Okamoto r-function theory of Py 
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2.1 The Jimbo-Miwa-Okamoto a-form of Py 



As formulated in |3C], the Okamoto r- function theory of the fifth Painleve 
equation Py is based on the Hamiltonian system {Q, P, t, K} 

(2.1) 

tK = Q{Q - ifp'^ - [{v2 - vi){Q - if - 2{vi + V2)Q{Q - 1) + tQ\p 

+ iv3 - Vl){v4 - Vi){Q - 1) 

where the parameters t^i , . . . , f 4 are constrained by 

(2.2) Vi+V2+V3 + Vi = 0. 



The relationship of (|2.l| ) to Py can be seen by eliminating P in the 
Hamilton equations 

^ ~ dP' dQ- 

One finds that Q satisfies the equation 
(2.4) 

y =(Tr + T)(y) -72/+ — 12 — [ay + -]+'y-+5 — 



with 
(2.5) 

11 1 

a = -^{v3-V4f, P = --{v2-vif, 7 = 2ui + 2t>2 - 1, '^ = "2" 

This is the general Py equation with ^ = (recall that the general 
Py equation with 5^0 can be reduced to the case with S = hy the 
mapping t \/—26t). Note that the first of the Hamilton equations ( |2.3D 
implies 

(2.6) 

tQ' = 2Q{Q - ifP -[{v2- vi){Q - 1)2 - 2{vi + V2)QiQ - 1) + ^q}, 

which is linear in P, so using this equation to eliminate P in shows 
that tK can be expressed as an explicit rational function of Q and Q' . 
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Of fundamental importance in the random matrix apphcation is the 
fact that tK (or more conveniently a variant of tK obtained by adding a 
term linear in t) satisfies a second order second degree equation. 



Proposition 2.1 With K specified by (|gJ[j, define the auxiliary 

Hamiltonian h by 



(2.7) 



h = tK + (^3 — vi){v4 — vi) — vit — 2v\. 



The auxiliary Hamiltonian h satisfies the differential equation 



(2.8) {th"f -(h-th' + 2{h'f)' + A\{{h: + vk) = 0. 



k=l 



Proof: Following |30|], we note from (|2.l|) , (|2.2D and the Hamilton 
equations (^^) that 



P 



h' = -QP - vi 

(2.9) th" = Q{Q^ - l)p2 _ [(^3 + v^- 2vi)Q^ -V2 + vi 

+ {v^ - vi){vi - vi)Q. 

We see from these equations that 

(2.10) th" = Q{h' + vz){h' + Vi) + P{h' + V2), 

while the first of these equations together with (^]^) and ( |2.7D shows 

(2.11) h-th' + 2{h'f = Q{h' + V3){h' + va) - P{h' + V2). 



Squaring both sides of ( p.lO ) and ( p.ll| ) and subtracting, and making 
further use of the first equation in (|2.9|) gives 



Actually the differential equation obtained by Jimbo and Miwa p^] is 



not precisely (2.8), but rather a variant obtained by writing 



h = a'^^'^ - Vjt - Ivj , 3 = 1,. 
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or equivalently 

(2.12) a(j') =tK + {v^ - vi){vi - vi) - {vi - Vj)t - 2{vl - v]). 

With this substitution ( p.^ ) coincides with the Jimbo-Miwa-Okamoto a 
form of Pv ( 1.37| ) with 



(2.13) {1^0,1^1, 1^2, l^s} = {Vl - Vj,V2 - Vj,V3 - Vj,V4 - Vj} 

(because ( |1.37| ) is symmetrical in {I'k} the ordering in the correspondence 
( p. 13 ) is arbitrary). 



The r-function is defined in terms of the Hamiltonian by 

(2.14) K=:^logT. 

at 



(2.15) ct(^') = t- log (e-^''^-^^)**^"^-"^)^"*""^)-'^"'-^' V) . 



It then follows from ( |2.12| ) that 

d_ 

'di 

Now, according to the Okamoto theory the r-function corresponding to 
some particular sequences of parameter values can be calculated recur- 
sively in an explicit determinant form. The above theory tells us that 
these T-functions can also be characterised as the solutions of a non- 
linear differential equation. The determinant solutions are a consequence 
of a special invariance property of the Py system which we will now sum- 
marise. 

2.2 Backlund transformations 



A Backlund transformation of ( |2.1| ) is a mapping, in fact a birational 
canonical transformation 

T{v;q,p,t,H) = {v;q,p,t,H), 

where "v^q^p^i^H are functions of q,p,t, H , such that the Py coupled 
system (p.3[) is satisfied in the variables {'v;q,p,i). Formally 



dp Adq- dH Adt = dp Adq - dH A dt. 
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In general, the significance of such a transformation is that it ahows an 
infinite family of solutions of the Py system to be obtained from one seed 
solution. From our perspective a key feature is the property 

(2.16) T-\H) = H 

which holds for some particular operators T which have a shift action on 
the parameter vector v = (vi, ^2, t'a, ^"4) whose components are referred 
to the standard basis. Explicitly, following [Q, we are interested in the 
transformation Tq which has the property ( p.l6| ) with the action on the 
parameters 



/ 1 1 1 3 

(2.17) To-V = [Vi- -,V2- -,V3- -,V4 + - 



1 1 1 3n 

Although the existence of Tq was established in ||3^ (termed the parallel 
transformation /(v)), as was the fact that To(t) = t, explicit formulas for 
Tq{P) and Tq(Q) were not presented. 

It is now realised |^, ^ that the required formulas are more readily 
forthcoming by developing the Py theory starting from the symplectic 
coordinates and Hamiltonian {q,p,t,H} specified by 

(2.18) 

tH := q{q - l)p{p + t) - {v2 - vi + V3 - V4)qp + {v2 -vi)p+ (f 1 - ^3)*^. 

The coordinates and Hamiltonians constituting the two charts and 
( ^.18| ) are related by the coordinate transformation [^] 

(g-l)(Q-l) = l 
(q - l)p + (Q - 1)P = V3-Vi 
(2.19) tK = tH + {V3 - vi){v2 - Vi) 



so in particular 1 — 1/q satisfies the Painleve V equation (2.4). 



The Hamiltonian ( |2.18 ) is associated with a set of four coupled first 



order equations in symmetric variables |2^ which admit an extended 
type A^r^^ affine Weyl group Wa = WaiA^^^) xi Z4 as Backlund transfor- 
mations. The generators of the group Wa = (sq, si, S2, S3, vr) obey the 
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algebraic relations 



(2.20) 



Si=l, SiSi±iSi = Si±iSiSi±i, SiSj = SjSi {j^i,i±l), 
vrsi = Sj+ivr, tt^ = 1, (i, j = 0, . . . , 3, S4 := sq)- 



The parameters of Py can also be characterised by the level sets of the 
root vectors ao, 01,02,03 (with oq + ai + 02 + = 6) of the A^^^ root 
system, and are dual to the vector parameter v in the following way 

ai{v) := V2- vi, a2{v) := vi - V3, a3{v) := V3 - V4 . 



In general by the convention of [27| the action of a group operation T on 
the roots is given in terms of its action on the components of v by 



(2.21) 



T{a) ■ V = a{T-^ ■ v) . 



Following [25 1 the action of the Backlund transformations on the root 



system a and coordinates q, p is given in Table 2.1. Note that with Tq 



defined by (|2.17D , 



To(a) = (ao - 1, 01,02, as + 1), 



so according to Table 2.1 we have 



(2.22) 



To = S3S2SITT, Tq ^ = vr ^SiS2S3. 



As is the case for the dynamical system associated with Piv there is a 
symmetric form for the system describing the Py transcendent [^], and 
is expressed in terms of the four variables 



(2.23) 



/o 
fi 

f2 



P + t 



Vt 

Viq 
P_ 

Vt 

h = Vt{i ■ 
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with the constraints /o + /2 = /i + /s = V^- Using these variables the 
Backlund transformations take the simple and generic forms of 

Si{aj) = aj - aijOi, vr(aj) = aj+i, 



(2.24) 



(y. ' 

Siifj) = fj + Uij-^, 7r(/j) = /j+i, 

Ji 



where the Cartan matrix A = {aij) and orientation matrix U = (uij) are 
defined by 



(2.25) A 




-1 

2 
-1 



-1 


-1 
2 



U 



1 



-1 







1 



-1 



Let us now compute the action of Tq ^ on the Hamiltonian tH. First 
we can easily verify from Table 2A that 

sqUH) = tH + ao— 1- ao("2 - 1) 

p + t 

s\{tH) = tH + ait + aias 

S2{tH) = tH — a2t + a2(«o — 1) 

S3(tH) = tH + aiQ3 
(2.26) TT{tH) =tH +{q-l)p- a2t. 

It follows from these formulas, Table |2.1| and ( 2.22| ) that 



(2.27) 



TQ^tH) =tH 



tH 



tH + qp. 



According to the relations in ( |2.18|) between tH and tK, q,p and Q,P, 
we therefore have 

T^HtK) =tK- Q{Q - 1)P + {vs - vi){Q - 1) 

(2.28) 

= tK =tK 

vh-+To-v cn— >Tq -a 

which is the result deduced indirectly in |^0|. Note from ( 2.19| ) that 



(2.29) 



1 



(v,-v,)-(T,-Hq)-l)T,\p) 
To\q)-l 
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So 

S3 

vr 



-QO 

ao + ai 
ao + as 



ai + ao 

— Ql 
"2 



02 

02 
02 + Ol 

—02 
02 + 03 

O3 



O3 
03 + Oo 

03 
03 + 02 
-03 

Oo 



01 



03 



f(g - 1) 



Table 2.1. 



Backlund transformations relevant for the 



Pv Hamiltonian ( 2.1 



where Tq ^{q) and Tq ^(p) can be deduced from Table 2.1 using ( 2.22| ) 



2.3 Toda lattice equation 

Introduce the sequence of Hamiltonians 

(2.30) K[n] := K 

VI— > (til — n/4,i)2 — n/4,t;3 — n/4,i)4+3n/4) 

and let r[n] denote the corresponding r-functions so that 

(2.31) K[n] = j^logT[n]. 

We have already remarked that a crucial feature of K[n] from the view- 
point of application to random matrix theory is the fact that it (or more 
precisely the quantity ( 2.12| )) satisfies the differential equation ( 1.37 ). An 



equally crucial feature for application to random matrix theory is the 
recurrence satisfied by {T[n]}. 



Proposition 2.2 p^ , \23^] The t -function sequence T[n\, corresponding 
to the parameter sequence {vi — n/4, V2 — n./4, ^3 — n/4, ^4 + 3n/4), obeys 
the Toda lattice equation 



(2.32) ^^logfM = + *:=4 

r"' n dt 
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where 
(2.33) 



f\n] := rV2g{^^4-^i+n)i 



T n 



Proof: Following |25| we consider the sequence of Hamiltonians 
H[n] associated with the given parameter sequence and the Hamiltonian 
( p.l8| ). Now, with f[n] denoting the corresponding r-function, it follows 
from ( ^l9|) that 



T n 



Since ( 2.32| ) is unchanged by the replacement f[n] i— > t"+^"f [n], it suffices 
to show that ( |2.32| ) is satisfied with f[n] replacing r[n] in ( |2.33| ). 
From the definitions 



(2.34) 
5 log 



f[n - l]f[n + 1] 



(TQ^{tH[n]) - tH[n]^ - (tH[n] - To{tH[n 



On the other hand, it follows from ( 2.271 ) that 



(2.35) (TQ^{tH[n]) - tH[n]^ - [tH[n] - To{tH[n])'^ 

= q[n]p[n]-Toiq[n])To{p[n]). 

With Tq^ specified from (|2.22| ), explicit formulas for To{q[n]) and ro(p[n]) 
can be deduced from Table |2.1|, and it can be verified using the Hamilton 



equations for H that the RHS of ( 2.35 ) is equal to 



(5 log {q[n]{q[n] - l)p[n] + {vi - V3)q[n] + (^4 - vi + n) 

= 6log^( tHWi] + (W4 -vi+ n)t] = 6 log ^6 log f e(''*-^i+")* 
dt.y / at \ 



T n 



dt y" ' ' "'"^ " dt 

Equating this with the LHS of (|2.34|) gives a formula equivalent to ( ^.32[) . 
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2.4 Classical solutions 



The Toda lattice equation ( p.32| ) is a second order recurrence, and so 
requires the values of f[0] and f [1] for the sequence members f[n], (n > 2) 



to be specified. It was shown by Okamoto |3C] that for special choices of 
the parameters, corresponding to the chamber walls in the underlying 
root lattice, the Py system admits a solution with t[0] = 1 and t[1] equal 
to a confluent hypergeometric function. 

Proposition 2.3 [3^ ] For the special choice of parameters 

(2.36) vi = va 

it is possible to choose r[0] = 1. Furthermore, the first member r[l] of 



the T-function sequence ( \2. 3j\ ) then satisfies the confluent hypergeometric 
equation 



(2.37) 



t(r[l])" + {V3-V2 + 1 + t)(T[l])' + (^;3 - ^^i)t[1] = 0. 



Proof: Write Q = Q[0],P = P[0]. In the case vi = V4 we see from 
(p.ip that it is possible to choose 



p = 0, K = 0, 



(2.38) 



the latter allowing us to take r[0] = 1. With ( 2.38| ) the equations (2.28) 
(|11) and d^) then give 



(2.39) 



tj^logT[l]{t) = {v3-vi){Q-l). 



The quantity Q must satisfy the Hamilton equation 
(2.40) 



tQ' 



dtK 



dP 



P-0 



-{V2 - vi){Q - If - {v^ - V2)Q{Q -l)-tQ 



where use has been made of (|2.2|). Substituting ( 2.39 ) in ( 2.40 ) gives 

O- ■ 



28 P.J. FORRESTER AND N.S. WITTE 

The two linearly independent solutions of ( 2.371 ) 

r^[l](t) =iFi{v3-Vi,V3-V2 + l;-t) 

(2.41) = ^ TT e-'''u'"-'"-\l - uY'-"' du 

^ ' T{v:i-vi)T{vi-V2 + l) Jo ^ ' 

and 

r^[l](t) = e-V(l +vi-V2,V3-V2 + l;t) 

(2.42) =— -/ e-'V-''^il + uy-'-''^-Uu 

(the superscripts "a" and "s" denote analytic and singular respectively, 
and refer to the neighbourhood of t = 0). 

As noted by Okamoto ||3^, it is a classical result that the solution of 
the Toda lattice equation in the case f [0] = 1 is given in terms of T[l](t) 
via the determinant formula 

(2.43) f [n] = det[5^'+'=r[l](t)],-fc=o,...,n-i. 

The formulas (|^ and (|2.42D substituted in (|]3|) with = vi and 
thus reading 

(2.44) f[n] = r'/2gnV[n] 

give the explicit form of f [l](t), and thus we have two distinct explicit 
determinant formulas for f\n]. 



2.5 Schlesinger Transformations 

We are now in a position to derive difference equations for the dynamical 
quantities of the Py system which are consequences of the Schlesinger 
transformations or shift operators corresponding to translations by the 
fundamental weights of the A^^^ lattice 

(2.45) To = S3S2Si7r, Ti = 7rS3S2Si, T2 = Sl7rS3S2, T3 = S2Sl7rS3. 



It is well known that these difference equations can be identified as dis- 
crete Painleve equations 1 34 1 satisfying integrable criteria analogous to the 
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continuous ones. Here we briefly demonstrate this and find a difl'erence 
equation for the Hamiltonians. 



Proposition 2.4 [341 The Schlesinger transformation of the Py system 
for the shift operator Tq~^ generating the parameter sequence 
(ao + n, 01,02,03 — n) with n £ Z, corresponds to the second order dif- 
ference equation of the discrete Painleve type, dPiy , namely 



(2.46) 



Vn 1 - 



n 



ynVn+l = t- 



Vn 

_ Xn + ^ai + ao -\- n 



n> 1 



where Xn = /oN/iN — \oli and yn = Vi/ fi[n]. With qp[n] = tH[n+l] — 
tH[n], the Hamiltonian times t, tH[n\, satisfies the third order difference 
equation 



(2.47) 

-tqp[n] =^{ao+n+qp[n])[tH[n\ + {a3-n)qp[n] - a2t] 

+ {a2-\-qp[n + l\)[tH[n] + {a3-n)qp[n] + (ao + ai+n)i]| 
-^^tH[n] + {a3-n)qp[n] + {a3-n-l){a2 + qp[n-\-l]) - a2t| 
x^{qp[n]-a3+n)[tH[n] + {l-ao-n)qp[n] + ait] 

+ {qp[n-l]-ai)[tH[n] + {l-ao-n)qp[n] - (02 + 03-''^)*]} 
H-|ti?[n] + {l — ao—n)qp[n] + {l — ao—n){qp[n—l]—ai) + aitj. 



Proof: The Schlesinger transformations for the shift operator Tq 
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that are relevant for the identification are 



(2.48) 



(2.49) 



(2.50) 



ToHfo) = f3-^ + 



ToHfi) = fo 



ao ao + "i+a2 
Jo 



h + 



ao+ai 



To 



Tt{h) = f2 + 



ao + ai 

h + -r 
Jo 

ao + ai 



QO 



Using ( 2.5C| ) in ( |2.48| ) along with the constraint to eliminate /s we can 
rewrite the latter equation as 



(2.51) /i[n] + — — = Vt-/o[n+l] + 



/oM 



Vt-/i[n+l]' 



This identity can then be employed in ( p. 49 ) to arrive at 



(2.52) /o[n]/i[n] + /o[n-l]/i[n-l] = ai + Vt/iN + 



n— as 



/iN, 



which is the first of the coupled equations ( ^.461) . The second equation of 
the set ( ^.461) follows immediately from ( |2.49| ) rewritten to read 



(2.53) 



1 



^ 1 /oN/iN+«o+^ 

/iN/i[»^+i] /oN/iN /oN/iN-ai 



The simplest way to derive the difference equation for tH\vi\ is to recast 
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the Schlesinger transformations in terms of canonical variables q,p, 

ao+ai+n 



(2.54) tq[n + l] = t + p[n] 



ao + n 

q[n\ + 



t + p[n] 



(an + n)t (n+l — as)t 

(2.55) p[n+l] = -tq[n] - ^ , + ^ ^ 



t + p[n] ao+ai+n 
p[n] H 



q[n] + 



(2.56) tq[n-l]= -p[n] + 



t + p[n] 
n — 03 1 — ao—n 



n 



(2.57) t + p[n-l]=tq[n] + 



02+03-"- 
q[n\ H 

as — n 

pin] + z r-r 

1 - q[n\ 

(02 + 03-"-)* 



03 — n 

p[n] + 



1 — q[n] 

Considering the first two equations ( p. 54 , 2^55] ) we find that 
(2.58) 

[ao + n + q[n]{t+p[n])][qp[n]{t+p[n]) - aip[n]+ a2t] 

qp[n+l\ + 02 = m r~FTi — 7 ^ 

p[n\[ai-q[n\{t+p[n\)\ + (oo+oi+njt 

and expressing this in terms of the Hamiltonian yields 
(2.59) 

[tq[n]+ao + n + qp[n]][tH[n] + {a3-n)qp[n]- a2t] 

gpln + 1 +02 = r- j -, r r--, z ^ — j— j z 

tH[n\ + {a3 — n)qp[n\ + (03 — n — l)ig[nj + (oq + 01 +n)t 

In an analogous way we find for the down-shifted product 
(2.60) 

[p[n\+a3-n-qp[n]][tH[n] + {l-ao-n)qp[n] + oit] 

(7pn — 1 — Ol = j— j 7 r j— j 7 r — j^-j 7 r-. 

til [nj + (1 — oq- n)(i'pH — {l — ao—n)p[n\ — {a2+as — n)t 
Now the up-shifted equation can be solved for q and the down-shifted one 
for p allowing their product qp[n] to be expressed in terms of just the 
Hamiltonian and other members of the qp sequence. The final result is 

O- ■ 
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3 Application to the finite LUE 
3.1 The case N = 1 

Comparing the case = 1 of the integrals ( |1.7| ) and ( |1.8| ) with those in 
( ^ ) and we see that 

£;i((0,5);a,/i)=C5"+^'+V^[l](s) 



^i((s,oo);a,/i) =Cs'^+^+V[l](s) 
RecalUng (|2l5| ), we thus have that 



^3 — ^1 — ^1+1 



— ^;2^ — a + 1 



f^log (t-^^i((0,t);a,^)) =a^«(t) 
t^log {t-^^mt,^)■a,^l))=a<'\t), 
where both a'^^^^ (t) and a^^^^ (t) satisfy the Jimbo-Miwa-Okamoto a form 



of Pv ( MTD with 

fo = 0, z^i = -n, f2 = a + 1, 1^3 = 1- 



Note that in the case ^ = this is consistent with (1.38). 



3.2 The general case 

Although it is not at all immediately obvious, the n x n determinant 
formed by substituting (g^) and in { ^A^j ) can be identified with 

the general n cases of the integrals ( |1.71 ) and ( |1.8D respectively. Consider 
first (|2.43| ) with initial value (|2.4lD . 

Proposition 3.1 Let f[n\ be specified by the determinant formula ( \2.4S ) 
with 

(3.1) f[l]it) = t^^''e\Fi{v3-vi,V3-V2 + l;-t). 

Then we have 
(3.2) 

n I 

f[n] oc n / dui e*"'np-^^(l - uiY^-"^^-^ J] (^^ " ^j)'- 
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It follows from this that 

(3.3) EN(.{t,ooy,a,fi) =Ct'^^+^'^^+^\[N]{t) 
and 

(3.4) log (t-^^'ENiit, oo); a, /x)) = C/^(t; a, fi) 

where Uj\f{t;a, fj,) satisfies the Jimbo-Miwa-Okamoto a form of the Py 
differential equation ( l-3'\ ) with 



(3.5) = 0, vi = -/i, V2 = N + a, = N, 

subject to the boundary condition 



j2 « + ^ 



(3.6) ^-(^'«'/^).~o"^ + ^ -^^^T^TT^*- 



Equivalently UN{t;a, ix) is equal to the auxiliary Hamiltonian { 2.11: ) with 
j = 1 and parameters (\3.3^). 



Proof: First observe that if {T[n]} satisfies the Toda equation ( ^.32 ) 



with t[0] = 1, then {i"''^r[n]} is also a solution which is given by the 



determinant formula ( ^.43|) with T[l](t) i— > t'^r[l](t). Choosing c = —1/2 
and substituting ( ^H] ) for t[l](t) we thus have 



j,k=0,...,n-l 



(3.7) t-"/2r[n] = det [6^+\Fi{vi-V2 + l,V3-V2 + l;t) 

where use has been made of the Kummer relation 

iFi(a,c; -t) = e~\Fi{c - a,c;t). 

Our first task is to use elementary row and column operations to eliminate 
the operator 5-'+'^ in (p. 7]). Starting with elementary column operations, 
in column k {k = n — 1, n — 2, . . . ,1 in this order) make use of the identity 



(3.8) 



6iFi{a,c;t) = a(iFi(a + l,c;t) - iFi{a,c;t)j 
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and add {vi — V2 + 1) times column k — 1. This gives 
t~'^/'^f[n] (X det\6^iFi{vi-V2 + l,V3-V2 + l;t) ... 

Si+k^\Fi{vi-V2 + 2,V3-V2 + l;t) 

Next, in column k (A; = n — 1, n — 2, . . . , 2 in this order) make use of the 



j = 0,...,n-l 
fe = l,...,n-l 



identity (3.8) again and add (f i — V2 + 2) times column A; — 1 to obtain 
(3.9) 

i""/^f[n] ocdet 5^ iFi{vi-V2 + l,V3-V2 + l;t) 

6^ iFi{vi-V2 + 2,V3-V2 + l;t) . . . 

P+''-^lFi{vi-V2 + 3,V3-V2 + l;t) . 



i = 0,...,n-l 
fc = 2,...,n-l 

Further use of (^^) in an analogous fashion gives 

(3.10) t-''/'^f[n](xdet\6^iFi{vi-V2 + l + k,V3-V2 + l;t)] 

At this stage we use elementary row operations to eliminate the oper- 
ation 6^ in ( 3.1C| ). In row j {j = n — 1, n — 2, . . . , 1 in this order) make 
use of the identity 

(3.11) ^iFi{a,c;t) = ^^^Fi{a, c + l;t) + iFi{a, c;t) 
at c 



and subtract row j — 1 to get 



det 



iFi (t>i -^2 + 1+^, "^3 -'^2 + 1; i) 

{v3 - vi - k)6^^^tiFi{vi-V2 + l + k,V3-V2 + 2;t) 



j = l,...,n-l 
fc=0,...,n-l 



Next, in row j (j = n — l,n — 2, . . . ,2 in this order) make use of the 
identity (|1|) again and subtract 2 times row j — 1 to obtain 



T n oc 



tdet 



iFi{vi-V2 + l + k,V3-V2 + l;t) 
{v3 - f 1 - k)tiFi{vi-V2 + l + k,V3-V2 + 2;t) 
iv3-vi - k)26^-'^t^iFi{vi-V2 + l + k,V3-V2+3;t) 



j = 2,...,n-l 
fc=0,...,n-l 
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^n(n- 1)/2 



{v3-Vi-k)jiFi{vi-V2 + l + k,V3-V2 + l+j;t) 



j,k=0,...,n—l 



thus eUminating entirely the operation S. 

We now substitute the integral representation for iFi deducible from 



( p.4lD to obtain 



j,k=0,...,n—l 



t 



n{n-l)/2 



n „i 



(3.12) 



X det 



j,k=0,...,n—l 

Regarding the determinant in the last line of this expression, we have 



det 



n _ 7, . , -, \n+j-k-l 



j,k=0,...,n~l 

n— 1 

= Jl (1 - Uj+iY det[u^+i]j- fc=o,., 

j=0 



,n— 1 ) 



where the equality follows by adding the k + 1-th column to the A;-th 
(A; = n — 1, . . . , Z) in a sequence of sweeps 1 = 1,... , n — 1. All factors other 
than the determinant in the multidimensional integral are symmetric in 
{uj} so we can symmetrise this term without changing the value of the 
integral (apart from a constant factor n! which is not relevant to the 
present discussion). Noting that 



n-l 



(3.13) Sym( - uj+iY det[uj^i]j^k=o,...,n-i) = — H (^i ~ ^^)^' 

j=0 ' j<k 



for some sign it, and substituting in ( |3.12| ) gives (|3.2|). 



r"/2f [nl oc t 



n(n-l)/2 



n 1 

n L dui e 

1=1-^^ 



j<k 
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Recalling (2.44), changing variables Uj i— > 1 — uj and comparing with the 
final integral in we thus have that 



T[N]{t) oc t-^'^+f'^^-^' ENiit, oo); a, /x) 
provided the parameters are given by ( |3.5D . The result (| 



now fol- 



lows by substituting this result in ( 2.15 ) with j = 1. For the boundary 
condition (^), we see from (|3.4|) and (|1.8|) that 



where Jn{o,, ij)[J2j=i ^j] denotes the integral ( |1.2C| ) with an additional 



factor of J2j=i ill the integrand. Noting that J2j=i can be written 
as a ratio of alternants allows the integral to computed using the method 
of orthogonal polynomials and leads to the result (^.6|). I 

Next we consider ( 2.43| ) with initial value ( 2.42 ). 

Proposition 3.2 Let f[n] be specified by the determinant formula ^-4^ 
with 

(3.14) f[l](t) = t^/^i;{vi-V2 + l,V3-V2 + l;t). 

Then we have 
(3.15) 



fin] oc 



n 



dui e 



Y[ {Uk - Uj) 

l<j<k<n 



It follows from this that 

(3.16) EN{{0,t);a,fi) = Ct^'^+^^^^+N' r [N] (t) 
and 

(3.17) log (t-^^^^((0, t); a, fi)) = VNit; a, /i) 



where Vi\i{t;a, fj,) is equal to the auxiliary Hamiltonian {2.1i) with j = 1 
and parameters ({3.^, and so satisfies the Jimbo-Miwa-Okamoto form of 
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Pv iHl-SV ) with parameters (3J.). The latter is to be solved subject to the 
boundary condition 

N{N + fi)a 



(3.18) F7v(i;a,/i) ~ 



Nt + N{a- fi) 
Proof: Analogous to (|3^) we have 



+ o(iA' 



(3.19) 



t 



-n/2 



T n 



det 



6^-^^^{vi-V2 + l,V3-V2 + l;t] 



' j ,k=0,...,n—l 

We now adopt the identical strategy as used in the proof of Proposition 
3.1, with the identities ( |3.8D and (3.11) replaced by 



6'ip{a, c; t) = a(^{a — c + l)ip{a + 1, c; t) — ip{a, c; t)j 



and 



di 



'ip{a, c; t) = tpi^a, c; t) — ■4'{a, c + l;t) 



respectively. The results (|3.15D and (|3.17| ) are then obtained by repeating 
the working which led to ( |3.2D and ( |3.4| ). 

For the boundary condition, we see from ( |3T7| ) and {^) that 



VN{t;a,fi) ■ 



-Nt + N{a - fj,) 



tIN{^J') 

where lN{lJ')[J2jLi ^j] denotes the first integral in ( |1.2C| ) with an extra 
factor of J2jLi ill the integrand. This integral can be computed by 
changing variables Xj ^ eXj in the definition of differentiating 
with respect to e, and setting e = 1. I 

In relation to (|1.15| ) we note from the fact that ( |2.37| ) is linear, with 
linearly independent solutions ( 2.41| ) and (2.42), that 

(3.20) 



r«[l](t) 



r(i + vi - V2) 

is proportional to the most general solution of ( p.37| ). Substituting ( ^.20 ) 
in ( p. 33 ) with n = 1 and V4 = vi, forming ( |2.43| ), and simplifying by 
noting that the considerations of the proof of Proposition 3^ (excluding 
the discussion of the boundary condition) remain valid independent of the 
value of ^, we obtain the following result. 
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Proposition 3.3 Let T[n] be specified by the determinant formula ( \2.4S ) 
with f [1] proportional to t^^"^ times /\3.20i ). Then 

(3.21) ^;v((0,5);a,/i;e) =Ct("+^)^+^V[iV](t) 
and 

(3.22) tj^ (t-^^^^((0, t); a, /i; O) = Ty^(t; a, ii) 

where W]\!{t;a, like UN{t;a,fi) from Proposition |g.i| and V]\i{t;a,id) 
from Proposition \3.^ is equal to the auxiliary Hamiltonian ( 2.1^ ) with j = 



1 and parameters (3.1), and so satisfies the Jimbo-Miwa-Okamoto form 
of PV ( 1-37 ) with parameters (3^). In the case /x = 0, with p{t) denoting 



the eigenvalue density of the L UE, the latter is to be solved subject to the 
boundary condition 

W7v(t;a,0) ^ T{N + a + l) 



(3.23) 



T{N)T{a + l)T{a + 2) 



e. 



3.3 A relationship between transcendents 



The evaluations ( |1.49D for Pmin(s; o) 
substituted into ( |1.11| ) imply that 



N^N+l 



and (|L40| ) for £^jv((0, s); a, 0) 



(3.24) VN{t;a,2) = -{a + l + 2N) +t + VN+i{t;a,{)) +t 



v^+i(^;«,o) 

VAr+i(t;o, 0)' 



In our previous study |17| we encountered an analogous identity relating 
some Piv transcendents. Like ( |3.24D , the Piv identity in [17| was discov- 
ered using a relation of the type (|1.11| ) . A subsequent independent deriva- 
tion using Backlund transformations was also found and presented. Like- 
wise the identity ( |3.24 ) can be derived from the formulas ( 2.26| ) and Table 
^!l| . An identical formula relates UN{t; a, 2) and U]\!^i{t; a, 0) which can be 
understood from the evaluation of pmax(s; a) and Eiy{{s, 00); a, 0). 

Proposition 3.4 Introduce the shift operators T2 and T3, which have 
the action on the parameters v specified by 



(3.25) 



T2 • V 



1 



1 



,V3 



,V4 



4'^ 4 
113^ 

{Vi - -,V2- -,V3 + 7,1^4-7 



4' 



4' 



4' 
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and consequently from ^.2^ and Table the representation 

(3.26) T2 = SITTS3S2, T3 = S2Sl7rS3. 

Then we have 

(3.27) TsT^^VN+iit;a,0) = FAr(t;a,2) 

and use of ({2. 2t ) and Table 2.1 on the LHS reduces this to ( 3.24 ). 

Proof: We take advantage of the arbitrariness in ( p.l3|) which with 
{ui} given by ( |1.38 ) with 1— > + 1 allows us to write 



(3.28) 



VN+i{t;a,0) = a 



(1) 



V2—vi=a+N+l 
V3—vi=0 
Vi—vi=N+l 



The actions ( 3.25| ) then imply 

r3r2-V7v+i(t;a,0) = a« 



V2—vi=a+N ) 
V3—vi=—2 
V4,—vi=N 



which recalling ( |1.41 ) and ( 2.13| ) is the equation ( |3.27| ). 

To evaluate the action of TsTg"^ on the LHS of ( |3.27| ), first note from 
( p. 12 ) and the final equation in ( ^.19| ) that for general parameters 



.(1) 



tH + (f3 - Vi){v2 - Vi). 



We remark that the simplifying feature of arranging the parameters so 
that a2 = ^3 — ^'i = in (|3.28 ) is that ( 2.18 ) takes the reduced form 



(3.29) 



tH 



Q2=0 



q{q - l){p + t) - (ai + 03)5 + ai 



The explicit formula ( 3.29| ) will become important later on. For now see 
seek the action of T3T2^ on a^^^ for general values of the parameters. 
First we note that the generators of the extended type A^^ affine Weyl 
group possess the algebraic properties ( |2.20D which can be used in the 
formula for T3T2~^ implied by (3.26) to show that 



TsT^ 



S2S3S1S0S1S3. 
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Now, it follows from ( 2.26| ) and Table 2.1 that 



(3.30) soSiS3{tH) =tH + ao ^ + ao(a2-l) + (ai+ao)*- 



Let us now consider separately the second term in (|3.30| ). Use of Table 
shows 

s^siiao-^) = (1-Q2)tg(g-1) 

+ t q{q - l){p + t) - ai{q - I) - a^q' 

Furthermore, we see from Table that in the special circumstance 02 = 
(i.e. ^3 — f 1 = 0), S2 acts like the identity on p and q so we have 

^ - 1) 

«2=0 q{q - l){p + t) - ai{q - I) - a^q 



S2S3Sl(ao 



p + t' 



(3.31) 



02=0 



where the final equality follows from ( 3.29 ) and ( 2.18| ). For the remaining 
terms in ( 3.30|) , we see from (IjH) and Table ^ that 



02=0 



S2S3S1 [tH + 00(02 - 1) + (oi +ao)*] 
The above results together with the simple formula 

2(l-ai) = 2(1-^2 + ^^1 



tH + t — oq. 



T3T2 H-ai"2^ 



02=0 



imply 

(3.32) T^T^^a^^^ 
= tH 



f3 — f 1=0 



+ t + I - {V2 - Vi) - {V4^ - Vi) + t- 
V3—v\=0 til 



^3— fi=0 



Substituting the values of V2 — vi and Vi — vi from ( 3.28 ) gives ( |3.24| ) 



3.4 Difference Equations 

We know the logarithmic derivatives of the r- functions Ui\i{t; a, ^) and 
V]\f{t]a, iJ.) satisfy the second order second degree differential equation 



( 1.37 ). Here we will utilise the Schlesinger transformation theory to show 
that they also satisfy third order difference equations in both a and /i 
variables. 
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Proposition 3.5 The logarithmic derivative U]\[{t; a, fi) satisfies a third 
order difference equation in the variable a 

(3.33) 

-t{U-U) = {{U - U+a+fi+l)[N{i2+t) + {a+l)U - aU] 

+ {a+l)t[V -U - N]} 

- [N{a+fi+l+t) + U - ia+l)iV - [/)} 

X ^^{U -U+a+fi)[Nfi + {N+a+fi+l)U -{N+a+fi)U] 
+ t[nN - fiU + {N+a+fi)U - {N+a)U]} 

- { - fj.{a+fi+t) + U- iN+a+n)(U - £)} 

where U := UN{t; a, fi),U := UN{t;a—l,ij),U := UN{t;a+l, fi), etc and 
the boundary conditions are expressed by UN{t',a,fj,) at three consecutive 
a-values for all N,t,fi. In addition U]\f{t; a, ^) satisfies a third order dif- 
ference equation in jj, 

(3.34) 

-t{U-Tf) = 

[iTf-U+2N+a+fi+l)[-N{N+a) + {N+fi+l)U - {N+fi)lf] 

- t[-N{N+a) + ifi+l)Tf - (iV+/i+l)lf + NU] } 

- { - N{2N+a+^l+l-t) + Tf- (A^+/i+l)(^ - [/)} 

X \^{Tf-tl+2N+a+fi)[-N{N+a) + {N+a+fi+l)U - {N+a+fi)Tf] 

- t[-N{N+a) -ntl + {N+a+ij)U - {N+a)lf] } 

- { - {N+a){2N+a+^^-t) + U - {N +a+ iJi){Tf -U)] 

where now U := UN{t;a, iJ, — l),Tf := U]s[{t;a, 



Proof: The difference equation generated by Tq^ ( 2.47 ) with the 
parameter identification V2 — vi = — //, v^ — vi = N + a,V4 — vi = N is not 
directly useful as this leads to a difference equation in both N and a (there 
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are no difference equations in the parameter alone - only combined 
ones with a or fi). However difference equations generated by the other 
shift operators can be simply found using this one with a permutation of 
the parameter identification. Thus for the difference equation in a one 
can use V2 — vi = —fi, V3 — vi = N,V4 — vi = N + a and the relation 



tH = U]sf{t; a, fi) + fiN which gives the result is (|3.33| ). The second result 
( |3.34| ) follows from the parameter identification V2 — vi = N, — vi = 
N + a,Vi-vi = -n and tH = UN{t; a, //) - N{N+a). I 

Both of these difference equations are of the third order and linear in 
the highest order difference and it may be possible to integrate these once 
and reduce them to second order equations, however we do not pursue 
this question here. Although we have stated the difference equations 
for UN{t;a, fi), it is clear that VNit]a,fj,) satisfies these as well although 
subject to different boundary conditions. 

3.5 ^Ar((0,s);a,^) for a G Z+ 

In a previous study the quantities E]\f{{0, s); a, 0) and E]\[{{0, s); a, 2) 
for a G Z+ were expressed in terms of a x a determinants. This was done 
using the method of orthogonal polynomials to simplify the corresponding 
multiple integrals. In fact we can easily express En{{0, s); a, /i) for a £ Z+, 
with fi general, as an a x a determinant using the methods of the present 
study. 



Proposition 3.6 For a G Z+ the function 

d_ 

'lit 

(3.35) 



at \ L i j,k=0,...,a-l'' 



dt 

Consequently, for a G Z>o, 



j,k=0,...,a-l 



r d^ 

(3.36) EnHO, s); a, /x) oc e"^^ det [—^^^^^(-s) 



j,fc=0,...,a-l 
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Proof: Choose 
(3.37) V3-V1 



-N, V3-V2 = fL 



in (3.1). Then according to (p. 431), for a £ Z> 



(3.38) 



-a/2 



T a 



det 



Sj+k^\Fi{-N,fi + l;-t) 



j,k=0,...,a-l 



Using the fact that 



iFi(-iV,/i + l;-t) cxL^(-i) 



where denotes the Laguerre polynomial, substituting ( 3.3^ ) in ( p. 441) 
and then substituting the resulting expression in (|2.15|) gives the first 
formula for cr^^^ in ( |3.35| ). In specifying the parameters in ( p. 15 ) we have 
made use of ( 3.37] ), the fact that f 4 — f 1 = a, as well as ( |2.2| ). The theory 
noted in the sentence containing (p.l3|) tells us that a^^^ satisfies ( 1.37 ) 
with 



2^0 = 0, Ui=V2-V3 = -H, V2=Vl-V3 = N, U3 = V4,-V3 = a + N. 

To obtain the second equality in (|3.35| ) we make use of the identity 

5 [L^(-i)e*] = [{N + l)L^+i(-t) - (Af + /i + l)L^(-t)] e* 

in conjunction with elementary column operations, proceeding in an anal- 
ogous fashion to the derivation of ( 3.101 ) using the identity (3^). This 
shows 



det 



oc det 



5^-(e*Li^^,(-t)) 



j,k=0,...,a—l L ' jy-rn,- ' ' i j.k=0,...,a—l 

The second equality in ( p.35| ) now follows by applying to this the general 
identities 



det 



6^{u{t)h{t)) 



det 



j,k=0,...,a-l 
j,k=0,...,a-l 



j,k=0,...,a-l 



dp 



jfkit) 



j,k=0,...,a-l 
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The function in the logarithm of the second equahty in ( |3.35| ) is of 
the form t~^^e~^^ times a polynomial in t. But for a G Z>o, we know 



from the second integral formula in (L7) that t~^^E{{0,t);a, fi) has the 
same structure. Since a^^^ also satisfies the same differential equation as 
VN{t;a,^) in ( |3.17| ), the formula ( 3.36 ) follows. I 



The determinant formula ( |3.38| ) can be written as an a-dimensional 
integral by using (^). Because ( ^.37 ) implies an exponent —{N + a) for 
the factors (1 — Uj) in the integrand, we must first modify the interval 
of integration. For this purpose it is convenient to first change variables 

I— > 1 — Uj. Then instead of the interval of integration [0, 1] we choose 
a simple, closed contour which starts at uj = 1 and encircles the origin. 
In particular, choosing this contour as the unit circle in the complex uj 
plane gives 



r"/2r[a] oc t 



a(a-l)/2 at 



Consequently we have the following generalization of (1.39). 
Proposition 3.7 For a G Z>o, 



CUEa 



(3.39) ^iv((0,s);a,/i) 
_jv, M,(0,0) 



CUEa 



We remark that the identity ( 3.3S| ) is itself a special case of a known 



more general integral identity The latter identity involves the PDFs 



N 



1=1 l<j<k<N 



Yl |e27r*x, _g2«x,|/3 (_i/2 < < 1/2) 

l<j<k<N 



APPLICATION OF THE r-FUNCTION THEORY 



45 



defining what we will term the ensembles L/3E and C/3E respectively. For 
/5 = 2, these PDFs were introduced in ( |1.4D and ( 1.32| ) as the LUE and 



CUE. For 13 = 1 and 4 these PDFs also have a random matrix inter- 
pretation. They correspond to the case of an orthogonal and symplectic 
symmetry respectively, and give rise to the matrix ensembles LOE, COE 
and LSE, CSE. Now define 

N 

4^)((0,.);a,M) = (nx[lU)(A^--)%^ 



(note that this reduces to (O) for (3 = 2). Then it follows from results 
in , derived using the theory of certain multi- variable hypergeometric 
functions based on Jack polynomials, that for a £ Z>o 

(3.40) E^^\iO,sy,a,,)=e--^ 7 

M^{2{^i+l)/f3-l,N) 

/il /c(4//3)eJ 

where M^^^ (a, h) denotes the integral ( |1.30| ) with exponent 2 in the prod- 
uct of differences replaced by j3. 

4 T-function theory of Pni and hard edge scahng 

In the Laguerre ensemble the eigenvalues are restricted to be positive. 
For N large, the spacing between the eigenvalues in the neighbourhood 
of the origin is of order By scaling the coordinates as in ( [L.54 ) 



the eigenvalue spacing is then of order 1 and well defined distributions 
result. In this limit the Py system degenerates to the Pni system so it is 
appropriate to revise the r-function theory of the latter. 



4.1 Okamoto r-function theory of P 



III 



Following [31 1, the Hamiltonian theory of Pni (actually the Pm/ system 



rather than the Pni ) can be formulated in terms of the Hamiltonian 
(4.1) tH = q^p^ - (g2 + viq - t)p V2)q- 
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Thus by substituting (0) in the Hamilton equations ( |2.3| ) and ehminating 
p we find that y{s) satisfies the Pm differential equation 

ds"^ y^ds'sdss y 
with q{t) = sy{s), t = and 

a = -4v2, P = 4(111 + 1), 7 = 4, 5 = -4. 



Note that with H specified by (|4.l| ), the first of the Hamilton equations 
(O) gives 



(4.2) tq' = 2q^p-{q^ +viq-t). 



Thus, by using this equation to eliminate p in (4J.), we see that tH can 
be expressed as an explicit rational function of q and q' . Analogous to 
Proposition it is straightforward to show that tH plus a certain linear 
function in t satisfies a second order second degree equation. 

Proposition 4.1 With H specified by define the auxiliary 

Hamiltonian 

(4.3) h = tH + ^vf-U 

The auxiliary Hamiltonian h satisfies the differential equation 

(4.4) {th"f + viV2h' - {A{h'f -l){h- th') - -^{vl + vl) = 0. 



Proof: Following [^], we note from (|4.1| ) and the Hamiltonian equa- 
tions ( |2^ ) that 

n = p 

(4.5) th" = 2(1 - p)pg + vip - -{vi + V2). 

Using the first equation to substitute for p in the second gives 

_th" -vih' +\v2 _ th" - vih' + \v2 

Ul-A{h'Y) ' (1 - 2/i') ' 
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while we can check from ( |4.3| ), ( [4.1| ) and the first equation in 

h-th' = {qp - ^vif - q[qp - ^{vi + V2)]. 
Substituting for qp and q, and simplifying, gives 
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that 



Of interest in the random matrix application is the variant of (4.4) 
satisfied by 



(4.6) 



(Tlll{t) := -{tH) 



vi, . t 
TV^l - "^2) + 7- 

t^t/i 4 4 



A straightforward calculation using the result of Proposition 4.1 shows 
that ajjj satisfies 



(4.7) 

- viV2{cf'iuf + (^'11 1 (Ml I - '^){crin - ta'jjj) 
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{vi-v2y 



0. 



Note that with the r-function defined in terms of the Hamiltonian ( [4.1| 
by ( 2.14| ), we have 



(4.8) 



at V 



4.2 Backlund transformations and Toda lattice equation 

For the Hamiltonian (4.1), Okamoto |31] has identified two Backlund 



transformations with the property ( p. 161) : 

(4.9) Ti-V = (V1 + 1,U2 + 1), T2-V = ivi + l,V2-l). 

The operators Ti and T2 can be constructed out of more fundamental 
operators so,si,S2 associated with the underlying B2 root lattice, whose 



action (following ||3l| and |^5[) on v,p,q is given in Table 4.2. According 
to Table t4.2| we have 



(4.10) 



Ti = S0S2S1S2, T2 = S2S0S2S1. 
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Vl 


V2 


P 


q 


t 


So 

Si 
S2 


-I-V2 
V2 
Vl 


— 1 — Vl 
Vl 
-V2 


f [qip-1) - livi-V2)\+l 
P 
1 — p 


t 

Q 

n 1 V2-V1 

-Q 


t 
t 

-t 



Table 4.2. Backlund transformations relevant to the 

Pin Hamiltonian (|0|). 



Analogous to the situation with the Hamiltonian ( 2.30| ) in the Py 
theory, introducing the sequence of Hamiltonians by 

(4.11) H[n] := H 

where H is given by ( |4.lD , the operator Ti can be used to establish a Toda 
lattice equation for the corresponding r- function sequence ( |2.31| ). 



Proposition 4.2 i31, 21] The r-funcUon sequence corresponding to the 
Hamiltonian sequence ( 4-11 ) obeys the Toda lattice equation 

f[n - l]f[n + 1] 



(4.12) 

where 
(4.13) 



6 log T[n] 



d 

S := t— 

dt 



f[n] :=r /V[n](t/4). 
Proof: Analogous to (|2.34D and (|2.35| ) we have 

(4.14) <^logltMi±Il = q[n]{l-p[n]) - T^\[n]{l - T^^p[n]). 
\n\ 

Making use of ( 4.10 ) we can use Table 4^ to explicitly compute T^^q\. 
and T^^p\p\. This shows 

q{n\(\-p\n\) - T{\\n\(\ - Tf ^M) 

= ~^(^^N/N - + wi)pN + 2(^1 + ^2))- 
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But according to ( |4.1| ) and the second equation in ( |4.5| ), this latter ex- 
pression is equal to 6log{-^tH). Substituting in ( 4.14| ) we deduce that 

d T[n - l]r[n + 1] 

—siogT[n] = c 2n 

at T^\n\ 



Substituting for r[n] according to ( 4.13| ) and taking C = 1 gives ( 4.12 ). 



4.3 Classical solutions 



Okamoto |31] has shown that for the special choice of parameters vi = 
—V2, the Piii system admits a solution with r[0] = 1, and allows t[1] to 
be evaluated as a Bessel function. 



Proposition 4.3 [31 j For the special choice of parameters 



(4.15) 



in (4jJ_) it is possible to choose r[0] = 1. The first member r[l] of the r- 
function sequence corresponding to ^.11 ), after the substitution t^t/A, 
satisfies the equation 



(4.16) t{f[l])" + {V, + l)(f[l])' - -f[l] = 0, f[l] := r[l](t/4). 

For Vl ^ Z, this equation has the two linearly independent solutions in 
terms of Bessel functions 



(4.17) 



f[l]=t-^^/2/±,^(^/t) 



(for Vl G I^^ and 1-^^ are proportional and (4jlj_) only provides one 
independent solution). 



Proof: Substituting (|4.15|) into the definition ( |4T| ) of H we see that 
it is possible to choose 



(4.18) 



p = 0, H = Q, 
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the latter allowing us to take r[0] = 1. To calculate q we use the Hamilton 
equation 

dtH 



(4.19) 



tq' 



dp 



p=0 



-{q' + viq-t). 



Now, it follows from ( |2.27] ) with H given by (U) and Tq^ by Ti that 



1 = 



Substituting this in ( 4.19 ) and changing variables t ^ t/A gives ( [4.16 ). 



For definiteness take the + sign in KW) and let vi = ly. This sub- 
stituted in ([4.13| ) with n = 1, and the corresponding formula for f[l] 
substituted in ( ^.431) shows that 



(4.20) 



^n{.-l)/2-[^] = det[5^'+^/,(^)],- fc=o,...,n- 



where use has also been made of the theory in the first sentence of the 
proof of Proposition |3.1| . The determinant in ( [4.20 ) can be written in a 
form which is independent of the operator 6. 



Proposition 4.4 We have 



(4.21) 



j,k=0,...,n~l 



n(n-l)/2 



det[Ij^k+u{Vi)]j^k=o,...,n-i- 



Proof: Let t = s"^ and note that 



1 



:= t— = -s— =: -Os 
dt 2 ds 2 

to conclude 

(4.22) det[S^+''UVi)],,k=o,...,n-i = 2-"("-^) det[<5^^/,(s)],, fc=o„ 



,,n-l- 



We now adopt a very similar strategy to that used in the proofs of Propo- 
sitions and |3.2|. Briefly, we first make use of the identity 



5sly{s) = Sly+l{s) + iylu{s) 
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together with elementary row operations to ehminate the operator 6^ in 
(|]22D, obt ainmff 



(4.23) det[5^+'=/.(s)],-fc=o,...,n-i = det[,5,^s^/.+,(s)]j- fc=o,...,n-i. 
Next, use is made of the identity 

(4.24) 5ssn,+,{s) = s^+i/,+,_i(s) - vsn,+j{s) 



to perform elementary column operations thus reducing the RHS of ( [4.23 ) 
to the form 



(4.25) 



det[s^I,+,-(s) . . . 5^is^+i/,+,_i(s) . . . ] .=0 



,ri-l . 
■ fc = l,...,n-l 



Next we apply elementary column operations, using the identity (|]2|) 
with J I— > J + 1, 1-^ — 2, to show that ( [4. 251) is equal to 



detis^L 



u+j 



.S] s- 



I j=0,...,n-l 
k — 2,...,n — l 



Continuing in this fashion, using the identity (4.24) with appropriate 
substitutions, reduces the RHS of ( [4.23| ) to 

det[s^^''I,y+j_k{s)]j,k=o,...,n-i- 
Substituting in (|^) gives (|42ll ). I 



Corollary 4.5 The function 



(4.26) aniit) = -tj log [e-'/H''"/' det[I,^k+AVi)U=o,. 



.,n-l 



satisfies the equation (\^. 1[ ) with parameters 
(4.27) {vi,V2) = {v + n,-v + n) 

and boundary condition 



^1/2 



(4.28) ani{t) . ~ , 7 " ^ " ( y " + 



t^oo 4 
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Proof: It follows from Proposition 0, ( [4.20| ) and ( [4. 131 ) that 

T[n]{t/A) oc t-^''^^det[Ij^k+uiVi)]j,k=o,. 



.,n-l 



The equation satisfied by (7jii{t) follows by substituting this, and the 
parameter values ( [4.271 ), i^i (i:^). To obtain the boundary condition, we 
make use of the well known Toeplitz determinant formula 



det 



1 

2^ 



j,k=0,...,n—l 



1 1 

n! (27r)^ 



1=1 



n/ deif{9i) n i^''" 



together with the integral representation 



Iniz) 



1 

2^ 



-in9+z cos 9 



l<j<k<n 



n G 



to rewrite the determinant in ( 4.26| ) as a multidimensional integral. 



(4.29) det[Ij^k+u{Vi)] 



j,k=0,...,n—l 



1 1 



n! (27r)" ^-^^^ 



n / dOie 



•/icos 6i—iu8i 



n 



l<j<k<n 



valid for i/ S Z. For large t the dominant contribution to the above integral 
comes from the neighbourhood of 9j = (j = 1, . . . ,n). Expanding to 
leading order about these points and changing variables shows 



t— >oo 



Substituting this in (|42^ ) gives (|^ 



The result of Corollary |4.5| is relevant to the hard edge scaling of 
En{{0, s); a, fi), which gives the quantity E^^"^{t; a, fi) defined by ( PTgoI) . 
Now, in an earlier study |]l6|, it was shown that for a £ Z>o 



(4.30) 



E'^^'^{t;a,0) = e-*/^det[/,_fc(^/t)],,fe=o,...,a-i 
E'^'''\t;a,2) oc e~*/^r" det[I,_fc+2( Vt)],-,fc=o,...,a-i 
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Thus, as already noted in [16| in the case of E^^"^{t; a, 0) (deduced from 
knowledge of ( |1.70 )), it follows from Corollary that E^^^'^{t;a, fi) for 
/i = and /X = 2 can be characterised as the solution of the equation (4.7) 
with parameters (f i, ^2) = (a + /x, a — /x). 

The results ( ^ ) were obtained by computing the limit of the RHS 
of (|3.36| ) using the asymptotic formula 

The same approach allows E^^'^'^{t; a, lu.) to be computed for general /x and 
a £ Z>o, giving the result 

(4.31) E''-'''{t;a,f,) oc e-*/^*-^"/^ det[/,-_,.+^(Vt)],,fc=o,...,a-i. 

It then follows from Corollary that the result ( |1.73| ) holds for general 
ji and a £ Z>o. For general a > — 1 the result (1.73) can be deduced from 
the first formula in ( p.. 40 ). This task will be undertaken in subsection 4.5. 
To conclude this subsection we make two remarks. The first is that 



substituting ( 4.29 ) in ( 4.31 ) shows that for /x G Z and a G 



^>o 



E^^-'-^it; a, /i) oc g-'/^t-Wa "Q f 

1=1'' ~ 



(4.32) 



cx e 



l<j<k<a 



This identity, relating an average in the infinite LUE scaled at the hard 
edge to an average in the CUE of dimension a, has previously been derived 
as a scaled limit of ( 3.39 ). With c := 2(/u + l)//3 and /x such that 
c G Z>o, it was also shown that the identity (3.40) in the hard edge limit 
reduces to 



(4.33) m---it-,a,,) = ce-^/^{;ty--^±y 



1=1 • 



l<j<k<a 



where 



^^-^r(l + 2//3)r(c + 2(j-l)//3) 



r(i + 2j/(3) 
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and E^l^')^^^'^^ like E^^^'^ specified by (1.69), is normalised so that it equals 
unity for t = 0. As a second remark, we note that a feature of the 
differential equation (4.4) is that it is unchanged by the mapping t —t, 
V2 ^ —V2- Using this, it follows from (^) and ( |4.6| ) that is also 

satisfied by 



(4.34) 



^niit) +1+^1^2/2 = -t;|log (e-^/H^^i^^-^-y^T{-t/4) ) 

where the equality follows from ( |4.8| ) . But by an appropriate modification 



d 



of Propositions 4.3 and KJ, a determinant formula for r(— 1/4) 



V2l-^-V2 

can be given for certain {vi,V2)- Performing these modifications and 



substituting the resulting formula in ( [4 .341 ) leads us to the conclusion 
that 

(4.35) -tj^ log (e-'/H"'/^ det[Jj^k+uiVi)]j,k=o,...,n-i 

satisfies ( [4. 7] ) with parameters 

{vi,V2) = {i^ + n,u - n). 



In the n = fi = 2 and v = a case, substituting ( [4.35 ) for (7{—t) in 
( |1.78| ) gives the well known jl^ expression for p^''^'^'^{s) in terms of Bessel 
functions. 



4.4 Schlesinger Transformations 

Integrable difference equations are found to arise from the Schlesinger 
transformations generated by Ti, T2 [|^, 33, 26] acting on q and p and we 
consider these here, along with the difference equations for the Hamilto- 
nians. 

Proposition 4.6 jg^/ The Schlesinger transformations of the Pni system 
for the operators Ti{T2) generating the parameter sequences {vi+n,V2+n) 
(respectively {vi+n, V2—n) ) acting on the transcendent q correspond to sec- 
ond order difference equations of the alternate discrete Painleve II, d-Pn , 
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equation (respectively its dual) 



l vi + V2 + 2 + 2n 1 vi + V2 + 2n _ _i q V2 + n 

^ ' 2 q[n]q[n + 1] + t ^ 2 q[n - l]q[n] + t ~ ^ t^ t ' 

1 V2- Vj - 2 - 2n 1 V2 - ri - 2n _ _i g fi +n 

^' 2 g[n]9[ri + 1] - 1 2 5[n- l]g[n] -t ^ t t ' 



Further, under the action ofTi, tH[n + 1] — tH[n] = q[n]{l — p[n]) and 
this satisfies the second order difference equation 



q{^-p)[n\[q{l-p)[n\ + ^(^i - ""2)} 

= t{q{l-p)[n + 1] + q{l-p)[n] + ^(t'l - ^^2)} 
(4.38) ^ i^q(^i_p^[n + l]+q{l-p)[n]+vi + l+n} 

X + q{i-p)[n - 1] + ^(vi - ^^2)} 

-j- {g(l-p)N + qil-p)[n — l]+vi + n|, 



while under the action 0/T2, tH[n+l] — tH[n\ = q[n]p[n\ and this satisfies 
the second order difference equation 



1 1 

QP[n][qp[n] - -{V1+V2)] =- t^qp[n + 1] + qp[n] - ^(^'i + ^^2)} 

■i- I qp[n + 1] + qp[n] — vi — 1 — n\ 
(4.39) ^ 

X |gp[n] + qp[n - 1] - -("Ui + V2)j 

■i- ^qp[n\ + qp[n — 1] — vi — n|. 



Proof: In the case of the Ti transformation the action on the canon- 
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ical variables in the forward and reverse directions are 



(4.40) q[n + 1] 



t 

— + 



Uvi + V2 + 2 + 2n)t 



q q[q{p - 1) - Un - V2)] + t 



(4.41) p[n + 1] = |[g(p -l)--{v,- V2)] + 1 



(4.42) q[n - 1] 



Uvi+V2 + 2n) 



P 



(4.43) p[n- 1] = 1- 



\{vi+V2 + 2n) 
P 

l{vi +V2 + 2n) 



(1 -P) - 2(^1 -^2) 



By combining (4.40) and ( [1.41| ) into the form 

t _ \{vi + V2 + 2 + 2n) 



(4.44) 



q\n + 1] + 



q[n\ 



p[n + 1] 



and ehminating p between this equation and (4.42) we have the result 
( [4.36| ). The difference equation for the Hamiltonian can be found by 
solving 



(4.45) q{l-p)[n + 1] + q{l-p)[n] + -{vi - V2) 



--{vi +V2 + 2 + 2n)— 1/ ^ — 



for q and 



(4.46) 



1 1 1 — 72 

q{l-p)[n - 1] + q{l-p)[n] + -{vi - V2) = -{vi + V2 + 2n) 



for p and then reforming (/(I — p). The result is ( [4.381 ). The corresponding 
results for the T2 sequence can be found in a similar manner. I 
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4.5 Hard edge scaling E^^^'^{t] a, fx) 

According to (|L69|) and ( |Oo|) we have for general a > — 1 



(4.47) 
where 
(4.48) 



E^'^'^'^is; a, fi) = exp / v{t; a, n) y 



{t;a,n) = hm ( Vat (t /4A^; a,fi) + ^n) . 



In this hmit the differential equation (1.37) characterizing the Py auxiliary 
Hamiltonian (2.12) and specifying Vn degenerates to the differential equa- 
tion (4.4) characterizing the Pni auxiliary Hamiltonian ( [4. 3D thus identi- 
fying V with this quantity. 

Proposition 4.7 The function v specified by ( 4-4^) satisfies the differ- 
ential equation 



(4.49) 



{tvy + a)'\v' f - v'{4v' + l){v- tv') - f^ifi + a)v' - ^ = 



42 



Consequently 
(4.50) 



v{t;a,fj,) 



crjiiit) + /x(/i + a)/2 



where am is specified by the Jimbo-Miwa-Okamoto a -form of the Painleve 
III equation ( 1.71 ) with parameters ( [ j . 7^ ) vi = a -\- fi,V2 = a — fi and is 
subject to the boundary condition ( 1 . 7l\ ) 



(4.51) 



i;(t;a,/i) ~ -\t + ]-at^/'^ -\a{a + 2^x). 

t^OO 4: 2 4 



Proof: Making the replacement a a — Nil in ( p37D with param- 
eters (|3.5| ), changing variables t 1— > t/AN, a{t/AN) h->- v{t), and equating 
terms of order N'^ (which is the leading order) on both sides gives (4.49). 
Substituting ( 4.50|) in ( 4.49| ) shows ani satisfies ( 1.71 ) with parameters 
as specified. The boundary condition is the same as ( |4.2g| ) with n = a, 
1/ = fi. I 
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Proposition 4.8 The function v{t; a, /x) satisfies a third order difference 
equation in the variable a (suppressing the additional dependencies) 

(4.52) ^a(a + l)t = [v{a+l) - v{a-l)+a+fi][v{a+2) - v{a)+a+fi+l] 

^t - av{a+l) + {a+l)v{a) 

and a third order difference equation in fi 

(4.53) -^t= [v{f,+ l) - v{fi)] [v{fi+l) - v{fi) + a] 

[vjn+l) - vip-l)+a+n][v{n+2) - v{n)+a+n+l] 
^ [v{fi+l) - v{n-l)+a][v{n+2) - v{n)+a] 

Proof: These two results follow from applying the hard edge scal- 
ing form ( [4.48| ) to the finite- difference equations ( p.33 ) and ( p. 34 ). 
Note that ([4.53| ) is precisely the result that would be inferred from the 
Pin difference equation generated by the T2 shift ( 4.39| ) with the corre- 
spondence 



(4.54) 



v{t; a, fi) = tH 



t/4. 



We have noted in ( |1.77 ) that a corollary of Proposition 4^ is the 
evaluation of Prmn{s', a) in terms of v{t; a, 2). Since, analogous to ( l.llj ), 

(4.55) ' ' 



P 



hard/ . ^ _ 
mm V^) - 



i?^-'i((0,s);a,0), 



the results ( |1.77| ) and ( [L.70D substituted into this formula imply an iden- 
tity between transcendents analogous to ( p^ . As is the case with ( p^ , 
this identity too can be independently verified. 



Proposition 4.9 With a{t) denoting the auxiliary Hamiltonian 
with parameters {vi,V2) = (a + 2,a — 2), and (TB{t) denoting the same 
quantity with parameters (^1,^2) = {a, a), the identity 



(4.56) 
holds. 



(j{t) = (JB{t) -l-t 
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Proof: The indirect derivation of this result has been sketched above. 
A direct derivation can be given by using the properties of the shift op- 
erator T2 in (|4.9| ) and ( 4.10 ). First, denote by tH[0] the Hamiltonian tH 



in (|4.6|) (and similarly define p[0], q[0]) so that 



aeit) = -{tH[0]) +t/A 



(4.57) 



It then follows from the definitions of T2 and a{s) that 



(4.58) 



a{t) = -{tTiH[0]) 



(a + 2) + t/4. 



On the other hand, from the definition ( [4.1| ) and the property ( 2.16 ) of 
T2 we have that 

tTlH[Q\ = tH[0] - q[0]p[0] - {T2q[0]){T2p[0]). 



But from (^lo| ) and Table U it follows that 



T2q[0] 
T2P[0] 



q[0] q[0]{q[0]p[0]-a)+t 



K#MO]-a) 



and thus 
(4.59) 



tTiH[0] = tH[0] -a-l + t 



where to obtain the final term use has been made of ( |4.1| ), ( [4. 61 ) and the 
analogue of the first equation in (^). Substituting in ( 4.5^ ) and recalling 
( ^1% gives (|]5|). I 



4.6 The spacing probability p2{0;s) 



The evaluation ( 4.47] ) has consequence regarding the nearest neighbour 
spacing distribution p2{0; s) for the scaled, infinite GUE in the bulk. This 
quantity has a special place in the Painleve transcendent evaluation of gap 
probabilities because it motivated the study of the probability £'2(6; s) 



60 



P.J. FORRESTER AND N.S. WITTE 



of no eigenvalues in an interval of length s in the scaled, infinite GUE 
through the relation 

cP 

^'2(0;^) = ^^2(0; s), 

and E2{0;s) in turn was the first quantity in random matrix theory to 
be characterised as the solution of a non-linear equation in the Painleve 
theory [^2|. More recently, building on the evaluation of E2{0;s) from 
[B^l, it has been shown that HH 



4.60 p2 0;s = ^ — ^exp / --^ dt 

s Jo t 

where a{s) is specified by the solution of the non-linear equation 

s\ay + 4{sa' - a){sa' -a + (a' f) - 4{a' f = 

(a close relative of ( 1.37| ) for a particular choice of the parameters) subject 
to the boundary condition a{s) ~ — (s^/Svr). Here we will provide a 
Painleve transcendent evaluation ofp2(0;s) distinct from ( [4.60| ). 
The starting point is the identity 

N 

(4.01) (n(>f-vx^,_,„,„.,>^,„, 

« £[A./2l((0, ^, 2)E[,«+,)/2l ((0, s'y,~,2), 

which is a variant of a formula given in , applicable to general matrix 
ensembles with a unitary symmetry and an even weight function. Now 
the LHS of ( [4. 611 ) multiplied by s"^ is proportional to the density function 
for a spacing between eigenvalues of length 2s symmetric about the origin. 
Recalling that the bulk scaling in the GUE requires 

ttXi 



A/ 



/2N 



it follows by replacing s by 7rs/v2iV, making use of the definition ( 1.69 ), 
and the property P2(0; s) ~ which follows from the boundary condi- 
tion for a{s) in ( 4.60[ ) (this fixes the proportionality constant) that 

(4.62) p2{0; 2s) = i^^hard((^^^2. 1 ^ 2)^t^-d((7rs)2; -1, 2). 
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According to ( [4.47D this specifies P2(0;s) in terms of Painleve III tran- 
scendents whereas ( [4.6C1| ) involves Painleve V transcendents. 

5 Concluding remarks 

The results of the paper have already been summarised in Section 1. 
Here we want to draw attention to a feature of this work (and our pre- 
vious work |l^]) which requires further study. This feature relates to the 
boundary conditions in the scaled limits, in particular the specification of 



E^^^'^{s; a, fi) by (1.73). One observes that the boundary condition (1.75) 



is even in /x, as is the differential equation ( |1.71| ) with parameters (|1.74 ). 



Thus according to this specification the solution itself must be even in ^u, 
but this contradicts the small t behaviour which for the formula ( [1.73 ) to 
be well defined must be 

(5.1) a{t) -i/i(/x + a) + 0{f), (6 > 0). 



The situation is well illustrated by the case a = 1, for which ( gl26|) 
shows 

(5.2) ^(i) = _tliog(e-*/W2j^(^) 
The small t expansion of I^{\/t) shows 

(5.3) -^'\Zo-r^^ + '^ + 4(iITT)'^ 

in agreement with ( ^.1]) , whereas the large t expansion of I^(\/t) shows 

where the a; are even functions of //. Thus the asymmetry between fi and 
— ;U can only be present in an exponentially small term for t ^ oo. This 
is clear from the exact expression 

sin vr/i 1 



(5.5) a{t;fi) - a{t;-fi) 



vr If,{Vi)I-^{Viy 
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Furthermore, although the boundary condition (|5.1| ) distinguishes the 
cases lb//, the differential equation ( 1.71| ) also requires that the leading (in 
general) non-analytic term also be specified for the solution to be uniquely 
specified (for example, with /i = this term is proportional to t°'~^^). 

A practical consideration of this discussion is that the boundary con- 
dition ( |1.75D does not uniquely determine the solution of ( |1.7lD , so that 
^hard^^. ^-j -g uniquely characterised. The same remark applies to 
the formula ( |1.60| ) for E'^°^^{s; ji). Indeed the inadequacy of the boundary 
condition ( |1.62D to uniquely determine the solution is already apparent 
in the our discussion of the formula (|l.63| ) . 

Acknowledgment. This work was supported by the Australian Re- 
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